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The auxiliary field method has been recently proposed as an elRcient technique to compute ana- 
lytical approximate solutions of eigenequations in quantum mechanics. We show that the auxiliary 
field method is completely equivalent to the envelope theory, which is another well-known proce- 
dure to analytically solve eigenequations, although relying on different principles a priori. This 
equivalence leads to a deeper understanding of both frameworks. 
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I. INTRODUCTION 



Finding analytical bound state solutions of the Schrodinger equation is a problem that deserves a constant interest 
in mathematical physics from the beginnings of quantum mechanics. Many different techniques have been developed: 
perturbation theory, WKB method, variational method, etc. We refer the reader to Ref. [l| for a complete review 
about analytical results in quantum mechanics. In this work, we would like to focus on two methods, both very 
efficient but less standard that the aforementioned techniques: the auxiliary field method (AFM) and the envelope 
theory (ET). 

The AFM has been first presented in Ref. [2] and extended in Ref. It allows to find an analytical approximate 
solution of the eigenequation H\^) = |^), with the standard Hamiltonian 

H = T(f) + V{r), (1) 

provided that the eigenequation 

Ha{v) \^a{v)) = [T{f ) + vP{r)\ \^a{v)) = €a{v) \^ a{v)) (2) 

is analytically solvable, v being a real parameter. The basic idea of the method consists in introducing an extra 
operator, called the auxiliary field, that transforms V{r) into a function of P{r) and of this auxiliary field. Analytical 
energy formulae can then be obtained by assuming that the auxiliary field is no longer an operator but a variational 
parameter. This is actually a mean field approximation with respect to the auxiliary field. The AFM has already 
been used to get accurate analytical energy formulae for radial power-law potentials of the form a + 6 r** 0, Q and 
for exponential potentials of the form —ar^ e~^'^ f3|. Notice also that, although these last three references deal with 
a nonrelativistic kinetic term of the form p'^ /{2m), analytical formulae can also be found in the case of a relativistic 
kinetic term of the form y'fp + rn^ [5] . 

The ET is nearly twenty years older than the AFM and has been proposed in Ref. [g']. The goal of this method 
is also to perform an analytical resolution of the eigenequation associated to the Hamiltonian ([T]). In the ET, the 
potential V{r) is replaced by a so-called tangential potential V*{r), depending on an extra parameter t and such 
that V{t -I- e) — V^{t -I- e) = O(e^) for any t and e <C 1. The tangential potential thus provides an approximation of 
the exact potential at each contact point r — t. One can say that it generates an envelope representation of V{r). 
Approximate analytical solutions for the eigenvalues associated with V{r) are eventually known provided that V*{r) 
is analytically solvable. The ET has currently been applied to a wide range of problems: funnel [3] and Yukawa 
potentials [1] , iV-body systems , relativistic Hamiltonians [13, EH , etc. 
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It has to be stressed that the AFM and the ET have been developed independently and are, a priori, two distinct 
frameworks from a conceptual point of view. In order to understand the peculiarities of these two methods, we present 
them in Sees. HI] and Uni Then, we show in Sec. IIVI that the AFM and the ET are equivalent, and we comment this 
equivalence. 



II. THE AUXILIARY FIELD METHOD 

We recall here the main points of the AFM and refer the reader to Refs. 0, 13] for more details. The problem is 
to find analytical approximate eigenenergies of the Hamiltonian ((T)). The AFM suggests the following procedure. We 
first assume that Ha = T(p^) + v P{r), where z/ is a real parameter, admits bound states with an analytical spectrum. 
The eigenequation ([2]) is thus analytically solvable, as well as the one associated with the Hamiltonian 

H{,y) = T{p^) + V{r) 

with V{r) = iyP{r)+V{I{i^))-i^P{I{v)), (3) 

the function I{x) = K^^{x) being the inverse function of K{x) defined by the relation (the prime denotes the 
derivative) 

V'(r) 

Kir) ^ (4) 

Let us remark that V{r) is of the form CiP(r) + C2 where Ci and C2 are constants. These numbers must be 
determined in order that V(r) approximates at best the potential V{r). Using the notation defined in Eq. ([2]), the 
eigenvalue E{v) of the Hamiltonian ([3]) is given by 

E{v)^eA{v)+V{I{v))-vP{I{v)). (5) 

The approximate eigenvalues and eigenstates of the Hamiltonian ([1]) are eventually given by E{vq) and \^ a{vo)) 
respectively, with achieving an extremum of the total energy ([5]), i.e. satisfying 

d,,E{v) = 0. (6) 

The value of vq depends on the quantum numbers of the state considered. Once z/q is known, the constants Ci and 
C2 can be computed. 

This procedure can be justified as follows. The arbitrary potential V{r) can be replaced by V{r,v) = v P{r) + 
V {!{£')) — V P {I{i>)), a function of the so-called auxiliary field i) and of an analytically solvable potential P{r). The 

exact elimination of this auxiliary field by S0V{r,i)) = gives i>o = K{r) and F(r, i>o) = V{r). Thus, the AFM 

would consequently lead to the exact results if the auxiliary field is seen as an operator i>o- As we assume it to be a 
variational parameter vq, the results are approximate but analytical. As intuitively expected, it can be checked that 
1^0 ~ (^'^(t'o)! i'o l^'yi(i^o)), with i>o = K{r) and given by ® [1]. The AFM can consequently be regarded as a 
"mean field approximation" with respect to a particular auxiliary field which is introduced to simplify the calculations. 

The main technical problem of the AFM is the determination of an analytical solution for the inversion of Eq. (2]) 
and for the determination of vq from Eq. Such a task can fortunately be achieved in many relevant cases t^^,,^,^]. 
It is of interest to mention here some general properties of the AFM [Sj] : 

• It is exact when V{r) — P{r). 

• The quality of the approximation can be analytically estimated. 

• The approximate eigenvalues are ruled by the correct scaling laws. 

• The accuracy of eigenvalues can be greatly improved, being still analytical, by a comparison with exact numerical 
results. 

• For a potential V{r) = vP(r) + (jv{r) where a is a small parameter and where av{r) can be treated in 
perturbation, {v{t)) = v{I{vq)) at the first order in a. 

• Let Pi{r) and P2{r) be two power-law potentials whose eigenenergies respectively read ei(A^i) and £2 (-^2), where 
Ni and A'^2 are the principal quantum numbers. Then, the eigenenergies of V{r) computed with both Pi{r) and 
P2{f) have exactly the same functional form, but depending on A'^i or on A'^2 respectively. 



III. THE ENVELOPE THEORY 



As the AFM, the ET Q is a method aiming to get approximate analytical energy formulae from an arbitrary 
Hamiltonian ([T]). We only present here its key features and refer the reader to Refs. [1, [13, [HI, [i2| for a detailed 
discussion about the basis and applications of ET. 

Let us set V(r) = v f{r) in Hamiltonian ([T]). Then the energy spectrum of this Hamiltonian is formally given by 
E — F{v), where the dependence on the usual quantum numbers n and I will be dropped for simplicity. The function 
F{v) is concave but not necessarily monotone. This allows to define a so-called kinetic potential k(S) by using the 
Legendre transformation (here the prime denotes the derivative with respect to v) 

k{S)^F\v) and S = F{v) - v F\v). (7) 

This transformation can be understood as follows. |^) being the eigenstates of Hamiltonian ([T]), one can define 
S = (*|r(p2) I*) and rewrite formally the energy spectrum as F{v) = S + v (*|/(r) |*) = S + vk{S). The 
transformation ([7]) follows from these relations. One is consequently led to the exact formula 

E ^ F(v) =mm\S + vk(S)]. (8) 
s>o 

What can now be done to go a step further in ET is to assume that V(r) — g{P{r)), where P{r) is a potential for 
which the solution of the eigenequation 

[T(p2) + y p(r)] \^a) = eA{v) \^a) (9) 

is analytically known. Then, 

{^a\T{p^)\^a) (10) 

can be analytically computed. It can moreover be shown that the kinetic potential corresponding to V{r), namely 
K{s), is given approximately by 

K{s)^g{kA{s)), (11) 

where ^^(s) is the kinetic potential associated to P{r). One then obtains an approximate form for the eigenenergies 
that reads [y] 

E £ = mm[s + g{kA{s))] . (12) 

s>0 

The variable s actually plays the role of a variational parameter. But, thanks to Eq. ©, the following equalities hold 

eAiv) = s + vkA{s), e'Aiv) = kA{s), (13) 
and another approximate energy formula coming from the rewriting of Eq. (|12p is 

£ = min [eA{v) - v e'A{v) + 5(4(«))] • (14) 

V 

This last formula is called the principal envelope formula in Refs. [l^, [lH ■ 

It is possible to understand Eq. (|14p as follows. If V{r) ~ g{P{r)), with g(x) a smooth function of x, then we can 
define the "tangential potential" V*{r) at the point r = < as 

W(r) = a{t) P{r) + g{P{t)) - a{t) P{t) 

with a{t)^^=g\Pit)). (15) 

Such a particular form is obtained by demanding that V^{r) and its derivative agree with V{r) and V'{r) at the point 
of contact r = t. If e ^ 1, one has indeed 

Vit + e)^v\t + e)^ '^P'itf g"{P{t)) + 0{e'). (16) 

The eigenenergies of Hamiltonian iJ* = T{p'^) + V^*(r), denoted by f (i), are given by 

£{t) - eA{a{t)) + g{P{t)) - a{t) P{t). (17) 



Let us now set 



t^a''^{v). (18) 
It can be computed from Eq. (|15p that a^^{v) = P^^{A{v)) with A{v) — g'~^{v), and Eq. ([T7]) becomes 



The final energy spectrum has to be minimized with respect to so we have 



and the physical energy reads 



£{v) =eA{v)+9iA{v))~vA{v). 
dvS{v)l^,^=Q^A{vo)^e'Avo) 
£{vq) = eA{vo) + g{£'^{vo)) - vq e^(wo), 



(19) 
(20) 
(21) 



that is nothing else than the principal envelope formula (|14p . 

We have just shown that ET can lead to analytical approximate energy formulae, namely Eqs. (fT2|) and (fH|) . which 
are both equivalent. Moreover, Eq. (|16p shows that £ is a lower (upper) bound of the exact energy if the function g is 
convex (concave), that is if g" > {g" < 0). The tangential potential indeed always underestimates (overestimates) 
the exact potential in this case. A clear advantage of ET is thus that it allows to know the variational or antivariational 
nature of the approximation that is performed. 



IV. EQUIVALENCE BETWEEN BOTH APPROACHES 

The similarity of the starting points of ET and the AFM is obvious: In both cases, a potential for which no 
analytical solution is known is "approximated" by an other potential for which analytical solutions exist. It suggests 
that a connection between both approaches should exist; and it will indeed be established in this section. Let us apply 
the AFM as described in Sec. |ll] with V{r) = g{P{r)). We find the following expression for the energy ([5]) 

E{u) = eA{v) + g {P{m)) - V P{IM), (22) 

the function I{x) — K~^{x) being computed from the relation ([4]). Remarkably, this AFM formula is equal to the 
ET one (fTO|) since I{x) = P~^{g'^^{x)). Consequently, the AFM and the ET lead to the same final energy formula 
(fT^ . The link between both approaches is given by 

1^ = a{t). (23) 

Moreover, with the point tq defined by the relation tq = the potential V{r) takes the form 

V{r) = K{ro) (F(r) - P(ro)) + l/(ro). (24) 

It is then easy to see that V(ro) = V{ro) and that V'{ro) ~ F'(ro). So, the potential V{r) is tangent to the potential 
V{r). An explicit example is presented in appendix. 

The function I{x) can be defined if the function K{x) can be inverted. In order to fulfill this condition, it is 
sufficient that K{x) is monotonic, that is to say that K'{x) has a constant sign. But, from the definitions above, we 
have K{x) = g'{P{x)), which imphes that K'{x) = g"{P{x)) P'{x). Since K{x) must be monotonic, the convexity of 
the function g is well defined if P{x) is also monotonic. This is the case if P{x) is a power-law potential, for instance. 
In these conditions, the convexity of the function g can also be used to determine the variational character of the 
AFM. Let us note that the existence of I{x) does not guarantee an analytical solution: The equation ^ must also 
be solved. 

Let us summarize our results. The auxiliary field v can be introduced as an operator in the Hamiltonian (HI), and 
leads to an equivalent formulation of this Hamiltonian. If one considers it as a variational parameter rather than an 
operator, as in the AFM, the results are approximate but can be analytical. We have shown in this section that the 
auxiliary field, when seen as a variational parameter, is nothing else than the function a{t) generating the tangential 
potential in ET. This shows that, although obtained in different ways, the AFM and the ET lead to the same results. 
In this way, some formulas about the power-law potentials obtained in Ref. [l3i | by the ET were rediscovered with the 
AFM in Ref. [2], but supplementary results are given in this last reference. 

Taking this equivalence into account, we can now better understand the meaning of the variational parameter v in 
the ET: Its optimal value will be close to (o'(P(r))) as clearly suggested by the AFM. Moreover, the properties of the 
AFM that have been proven in Refs. also hold for ET. Finally, we can now have an a priori knowledge of the 

(anti)variational nature of the AFM energy formulae provided that we express V{r) as g{P{r)) and compute whether 
g is convex or concave. 



Acknowledgments 



C. S. and F. B. would thank the F.R.S.-FNRS for financial support. The authors thank the anonymous referee for 
useful suggestions. 



APPENDIX A: THE ANHARMONIC OSCILLATOR 



In order to illustrate both methods described above, they will be applied to the case of the 3D anharmonic potential 
V{x) = 3 .T^ + 8y/P X for the following dimensionless Hamiltonian 

H + 3x^ + 8y^x, (Al) 

where q and x are conjugate variables. This example is chosen because technical details can be found in Ref. [3]. In 
this paper, it is shown that, using the AFM with the potential P{x) — x^, the approximate eigenenergy of H, with 
quantum numbers n and is given by 



eapp(/3;n,0 = 2/3y (^Gi(r) + ^^j with Y = l^j j , (A2) 

where G- is a complicated function explicitly given in Ref. [Sj and where N — 2n + I + 3/2. 

Following the discussion above, the ET would give exactly the same form for eapp. But, in addition, with this 
method, it is now possible to show that the function g{y) = 3y + 8y/Py is concave {g"{P{x)) — —2y/P\x\^'^ < 0) so 
that Eapp is an upper bound of the exact results. This has been checked numerically for various values of (3. The 
envelope potential of V{x) is then given by 

V{x)^ (^3 + ^^x^ + 4^xo, (A3) 
where the number a;o depends on the state considered. More precisely, at the minimization point, 

Xo = with = -p—^. (A4) 

It is easy to see that V{xo) = V{xo) and that V'{xo) — V'{xo). Moreover, the relation 

V{x)~V{x) = ^{x-xof >0 (A5) 

Xo 

guarantees that eapp is always an upper bound. 
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